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ABSTRACT 

Using large deviations results that characterize the amount 
of information per node on a two-dimensional (2-D) lattice, 
asymptotic behavior of a sensor network deployed over a 
correlated random field for statistical inference is investi- 
gated. Under a 2-D hidden Gauss-Markov random field 
model with symmetric first order conditional autoregres- 
sion, the behavior of the total information [nats] and energy 
efficiency [nats/J] defined as the ratio of total gathered in- 
formation to the required energy is obtained as the coverage 
area, node density and energy vary. 

1. INTRODUCTION 

In this paper, we investigate the fundamental behavior of a 
flat multi-hop ad hoc sensor network deployed over a cor- 
related two-dimensional (2-D) random field for statistical 
inference. In particular, we examine the amount of informa- 
tion obtainable from a sensor network distributed over a 2-D 
Gauss-Markov random field (GMRF) and related trade-offs 
in various asymptotic settings. We consider the Kullback- 
Leibler information (KLI) and mutual information (MI) H] 
as our information measures. Our approach to calculating 
the total obtainable information is based on the large devi- 
ations principle. That is, for large networks the total infor- 
mation is approximately given by the product of the number 
of sensors and the asymptotic per-sensor information. How- 
ever, a closed-form expression for the asymptotic per-sensor 
information (or asymptotic information rate in 2-D) is not 
available for general 2-D signals. To address this problem, 
we adopt the conditional autoregression (CAR) model and 
corresponding correlation model for the signal, and derive a 
closed-form expression for the asymptotic information rate 
in 2-D. We do so by exploiting the spectral structure of the 
CAR signal and the relationship between the eigenvalues of 
the block circulant approximation to a block Toeplitz matrix 
describing the 2-D correlation structure. Based on the de- 
rived expressions for asymptotic information rate and their 

^ Y. Sung and H. Yu are with the Dept. of Electrical Engineering, KAIST, Daejeon 
305-701, South Korea. Email:ysung@ee.kaist.ac.kr and hjyu@stein.kaist.ac.kr. H. V. 
Poor is with the Dept. of Electrical Engineering, Princeton University, Princeton, NJ 
08544. Email: poor@princeton.edu. The work of Y. Sung was supported in part by 
Brain Korea 21 Project, the School of Information Technology, KAIST. The work 
of H. V. Poor was supported in part by the U. S. National Science Foundation under 
Grants ANI-03-38807 and CNS4)6-25637. 



properties, we investigate the behavior of sensor networks 
deployed over correlated random fields for statistical infer- 
ence. 

1.1. Related Work 

Large deviations analysis of Gauss-Markov processes in Gaus- 
sian noise has been considered previously. (See 111 and 
references therein.) However, most work in this area con- 
siders only one-dimensional (1-D) signals or time series. 
A closed-form expression for the asymptotic KLI rate was 
obtained and its properties were investigated for 1-D hid- 
den Gauss-Markov random processes Q. Large deviations 
analyses were used to examine the issues of optimal sensor 
density and optimal sampling in a 1-D signal model in jSl 
and iHj. For a 2-D setting, an error exponent was obtained 
for the detection of 2-D GMRFs in ||5l, where the sensors 
are located randomly and the Markov graph is based on the 
nearest neighbor dependency enabling a loop-free graph. In 
this work, however, measurement noise was not considered. 
Our work here focuses on the analysis of the fundamental 
behavior of 2-D sensor networks deployed for statistical in- 
ference via new large deviations results for 2-D hidden GM- 
RFs, which enable us to investigate the impact of field corre- 
lation and measurement signal-to-noise ratio (SNR) on the 
information. 

2. BACKGROUND AND SIGNAL MODEL 

To simplify the problem and gain insights into behavior in 
2-D, we assume that sensors are located on a 2-D lattice 
X„ = [0 : 1 : n — 1]^, as shown in Fig. [T] We assume 
that the signal samples of sensors form a (discrete-index) 2- 
D GMRF and that each sensor has Gaussian measurement 
noise. The (noisy) measurement of Sensor ij on the 2-D 
lattice In is given by 

y„- = x„- + ij € J„, (1) 

where {Wij} represents independent and identically dis- 
tributed (i.i.d.) A/^(0, cr^) noise with a known variance tr^, 
and {Xij} is a GMRF on the 2-D lattice independent of 
the measurement noise {Wij}. Thus, the observation sam- 
ples form a 2-D hidden GMRF. In the following, we briefly 
introduce the results on GMRFs relevant to further develop- 
ment. 



Definition 1 (GMRF lO) A random vector X = (Xi, X2, 
• • • ^Xn) G M" is a Gauss-Markov random field with re- 
spect to (w.r.t.) a labelled graph Q ~ ii^,^) with mean fi 
and precision matrix Q > 0, if its probability density fiinc- 
tion is given by 



p(X) = (277)- 



/2|Q|l/2 



cxp 



1 



-(x-m)^Q(x-m) 



(2) 



and Qirn 7^ •i=^ {/, m} G £ for all I 7^ m. Here, v is 
the set of all nodes {1, 2, • • • , ??} and £ is the set of edges 
connecting pairs of nodes, which represent the conditional 
dependence structure. 




Fig. 1. Sensors on a 2-D Lattice X„: Hidden Markov Struc- 
ture 

The 2-D indexing scheme ij in ([T]) can be appropriately 
converted to an 1 -D scheme to apply Definition [T] From 
here on, we use the 2-D indexing scheme for convenience. 

Definition 2 (Stationarity) A GMRF {Xy } on a 2-D dou- 
bly infinite lattice X^o is said to be stationary if the mean 
vector is constant and Gov{Xij , Xi>ji) = E{(Xy— E{Xy }) 
{Xi>j> — E{Xi/j/})} = c(i — — j') for some function 
<;■). 

For a 2-D stationary GMRF {Xij }, the covariance {7^ } 
is defined as 7^ = E{Xj'jvXj'+i.j/+j} = E{XooXy}, 
which does not depend on i' or j' due to the stationarity. 
The spectral density function of a zero-mean and stationary 
GMRF on Too with covariance 7^ is defined as 

f{uji,uJ2) = ^ ^ 7yexp(-t(icji + jcJ2)), (3) 



wheret = V^and {u;i,lu2) G (— 7r,7r]^. Note that this is a 
2-D extension of the conventional 1-D discrete-time Fourier 
transform (DTFT). 

Definition 3 (The Conditional Autoregression ) A GMRF 
{Xij} is said to be a conditional autoregression ( CAR) if it 
is specified using a set of full conditional normal distribu- 
tions with mean and precision: 



E{Xy|X_y} 



1 



Prec{Xij\K.ij} = 000 >0, (5) 
where ^-ij denotes the set of all variables except Xij. 



It is shown that the GMRF defined by the CAR model Q - 
Q is a zero-mean stationary Gaussian process on Too with 
the power spectral density Q 



(6) 



if |{6»ij 5^ 0}| < oo, =e_.i,_j, 6»oo>0, (7) 
{6»ij}issothat/(a;i,(.U2) > 0, \/{loi,W2) S {-n,iTf. (8) 

Henceforth, we assume that the 2-D stochastic signal {Xij } 
in ([T]i is given by a stationary GMRF defined by the CAR 
model dUi - © and dill - (l8]l. 

3. ASYMPTOTIC INFORMATION RATES AND 
THEIR PROPERTIES 

In this section, we derive a closed-form expression for the 
asymptotic KLI rate and MI rate in the model ([T]i, defined 

as 

OC = lim log— {{Yij,ij € I„}) a.s. under po, and 

ri^oo |i,„| pi 

I = lim -^I{{Xij,ij e Tn}; {Yij,ij e J„}), 

respectively. For the MI, the signal model ([T]i is directly ap- 
plicable, whereas for the KLI the probability density func- 
tions of the null (noise-only) and alternative (signal-plus- 
noise) distributions are given by 



PoiYij) : Yij=Wij, ijeTn, 



^X^J+W^,, tjeTn. (9) 



The following closed-form expressions for the asymptotic 
information rates in the spectral domain have been obtained 
in Q by exploiting the spectral structure of the CAR sig- 
nal and the relationship between the eigenvalues of block 
circulant and block Toeplitz matrices representing 2-D cor- 
relation structure. 

Theorem 1 For the model (O with the signal given by (|4| - 
(O, assuming that conditions (l7| - (O hold, the asymptotic 
KLI rate is given by 



3C = 



log 



(10) 



dujid(jj2 



2 + 4:Tr^f{ui,uJ2) 2, 
= J-^ I I D(J\r(0,Sl(u:i,ui2))\W(0,Sl(aii,uJ-2)) du:xdu:2 

where D{-\\-) denotes the Kullback-Leibler divergence. 
Proof: In 0. 



As a by-product of the proof of the above theorem, we 
have the asymptotic MI rate given by 



Corollary 1 The asymptotic KLI and MI rates for the SF- 
CAR 2D signal model are given by 



1 

4^ 



47r2/(tJi,CJ2] 



dLUiduj2- 



(11) 

Theorem[T]is a 2-D extension of the asymptotic KLI rate of 
1-D hidden Gauss-Markov model obtained in and the 
asymptotic KLI rate ( fTOb can be explained using a frequency 
binning argument. Specifically, for each 2-D frequency bin 
dujiduj2, the spectra are flat, i.e., the signals are independent 
and Stein's lemma can be applied for the bin. The overall 
KLI is the sum of contributions from each segment. 

3.1. Symmetric First Order Conditional Autoregression 

To investigate the properties of the asymptotic KLI and MI 
rates as functions of field correlation and SNR, we further 
consider the symmetric first order conditional autoregres- 
sion (SFCAR), defined by the conditions 



E{X,,|X_,,} 



A 



+ x.. 



PKc{Xij\X-,.j} = K>0, 



where < A < -1. (This is a sufficient condition to sat- 
isfy (|7]i - ([8]l.) Here, ^oo = k and 0i,o = 9^1^ = 0o,i = 
(?o,-i = —A. In the SFCAR model, the correlation is sym- 
metric for each set of four neighboring sensor nodes. The 
SFCAR model is a simple but meaningful extension of the 
1-D autoregression (AR) model which has the conditional 
causal dependency only on the previous sample. Here in 
the 2-D case we have conditional dependence on four neigh- 
boring nodes in the four (planar) directions, capturing 2-D 
correlation structure. The spectrum of the SFCAR signal is 
given by 



1 



47r2K(l — 2^cosa;i — 2(;^coscj2) ' 



where the edge dependence factor C is defined as 

C = -, o<c<i/4. 

K 



(12) 



(13) 



Here, C = corresponds to the i.i.d. case whereas ( — 1/A 
corresponds to the perfectly correlated case. Therefore, the 
correlation strength can be captured in this single quantity 
C for SFCAR signals. The power of the SFCAR is obtained 
using the inverse Fourier transform via the relationship (O, 
and is given by Ps = 700 = (O < C < i). where 

A'( ) is the complete elliptic integral of the first kind ||9l- 
The SNR is given by SNR = ^ = Using ([T0| and 

the SNR, we obtain the asymptotic KLI and MI rates for the 
SFCAR signal, given in the following corollary to Theorem 
1, also from Q. 



1x2 ./-TT ./-X V 2 y (2/7r)it(4<)(l - 



SNR 



2C cos 1^1 - 2C cos 1^2)/ 
(14) 



(2/x)iC(4C)(l-2C ooBuJi-2( cos 1^2) 

and 

's = /" - 1°S ( 1 + — 

4^2 J-ir J-TT 2 V (2 

respectively. 



(2/x)if(4C)(l - 2C costal - 2C cos 1^2) 



Note that the SNR and correlation are separated in ( fT4l i- 
(flST l. which enables us to investigate the effects of each term 
separately. 

3.2. Properties of the asymptotic KLI and MI rates {Xg 
and Is) 

First, it is readily seen from Corollary [T] that %s and Is 
are continuously differentiable functions of the edge de- 
pendence factor (0 < C 1^ 1/4) for a given SNR since 
f : X ^ K{x) is a continuously differentiable function 
for < a; < 1 flOl . The values of JCs at the extreme correla- 
tions are given by noting that A'(0) = ^ and K{1) — 00. 
Therefore, in the i.i.d. case (^ = 0), the corollary reduces 
to Stein's lemma as expected, and OCs is given by 



; log(l+SNR)H 



2(1 + SNR) 2 



^ r>(A^(0, l)\\^f{0., 1+SNR)). 



In the i.i.d. case, the asymptotic MI rate is given by the well 
known formula, /s I ij^o = ^log(l+SNR). For the perfectly 
correlated case (C, ~ 1/4), on the other hand, 3Cs = and 
Is = 0. (In this case as well as in the i.i.d. case, the two- 
dimensionality is irrelevant.) The limiting behavior of the 
asymptotic information rates is given by Taylor's theorem. 
Due to the continuous differentiability, we have 

XsiC) = ci .(1/4-0+0(11/4- CI), (16) 
/,(C) = 4 .(1/4-0+0(11/4- CI), (17) 

for some constants ci and c[, as C ^ 1/4. Similarly, we 
also have the linear limiting behavior for OCs and Ig in a 
neighborhood of C = with non-zero limit values, as C ^ 
0. That is, 

%siO = 3C,(0) + C2C + o(0, (18) 
/,(0 = /.(0)+c'2C + o(0, (19) 

for some C2 and c'a, as C — * 0. For intermediate values of 
correlation, it is seen that at high SNR OCs is monotonically 
decreasing as C increases. At low SNR, on the other hand, 
correlation is beneficial to the performance. 

With regard to %s and Is as functions of SNR, the be- 
havior of 3Cs is given by the following theorem. 



Theorem 2 The asymptotic KLI rate OCgfor the hidden SF- 
CAR model is continuous and monotonically increasing as 
SNR increases for a given edge dependence factor < 
^ < 1/4. Moreover, OCs increases linearly with respect to 
^ log SNR as SNR — > oo. As SNR decreases to zero, on 
the other hand, OCs converges to zero with the convergence 
rate %s{SNR) = C3 • SNR'^ + o{SNR^) for some constant 
C3 as SNR 0. The asymptotic MI rate Is has similar 
properties as a function of SNR, i.e., it is a continuous and 
monotonically-increasing function of SNR. At high SNR, it 
increases with rate ^ log SNR, whereas it decreases to zero 
with rate of convergence Is{SNR) — C3 • SNR + o{SNR) for 
some constant C3 as SNR —>■ 0. 

Proof: In ©. 

Note that the limiting behavior as SNR ^ is different 
for OCg and Ig', OCs decays to zero quadratically while Is 
diminishes linearly. 



The correlation function ( I2TI ) can be regarded as the repre- 
sentative correlation in 2-D, similar to the exponential cor- 
relation function e~^'^" in 1-D. Both functions decrease 
monotonically w.r.t. (i„. However, the 2-D correlation func- 
tion is flat at dn = ifTTl . Further, we have a continuous 
and differentiable mapping g : p C, from the edge cor- 
relation factor p to the edge dependence factor given by 

El 



p 



(2/^)if (40 - 1 



(23) 



4(2/7r)Ci^(4C) 

which maps zero and one to zero and 1/4, respectively. Thus, 
we have C — 9{f{dn)), and for given physical parameters 
(with a sUght abuse of notation), 

3C,(SNR, - 3C.(SNR, = 3C,(SNR, d„). 

(And, similarly for Is.) We will use the arguments SNR and 
( for 3Cs and Is properly if necessary. 



4. SCALING LAWS IN AD HOC SENSOR 
NETWORKS OVER CORRELATED RANDOM 
FIELD 

Based on the results in the previous sections, we are now 
ready to answer some fundamental questions in the design 
of sensor networks for statistical inference about the under- 
lying stochastic field. 

4.1. Physical correlation model 

The actual physical correlation for the SFCAR model is 
given by solving the corresponding continuous-index 2-D 
stochastic differential equation (the stochastic Laplace equa- 
tionf] El 



oxj \dyj 



X{x,y) = u{x,y), 



(20) 



where u{x, y) is the 2-D white zero-mean Gaussian pertur- 
bation and a > is the diffusion rate. By solving the SDE, 
the edge correlation factor p is given, as a function of the 
sensor spacing (i„, by ifTTI 

P = — = — = .f{dn) = adnKi{adn), (21) 
700 700 

where ) is the modified Bessel function of the second 
kind whose asymptotic behavior is given by 



Ki{x) 



as a; 00, 
1/x as x 0. 



(22) 



'Note that the solution of )20t is circularly symmetric, i.e., it depends 
only on r = ^x'^ + y'^, and samples of the solution y) of j20t on 
lattice In do not necessarily form a discrete-index SFCAR GMRF. How- 
ever, {2Q\ is still the continuous-index counterpart of the SFCAR model, 
and we use its correlation function for the SFCAR model. 



4.2. Asymptotic behavior 

In the following, we summarize the assumptions for the pla- 
nar ad hoc sensor network that we consider 

(A.l) sensors are located on the grid X„ = [0:1: n— 1]^ 
with spacing dn, as shown in Fig. [T] and a fusion 
center is located at the center ([n/2j , [n/2j ). 

(A. 2) The observations {Ky } at sensor nodes form a 2-D 
hidden (discrete-index) SFCAR Gauss-Markov ran- 
dom field on the lattice for each dn > 0, and the edge 
dependence factor is given by i2l[ and ( |23] |. 

(A. 3) The fusion center gathers the measurement from all 
nodes using the minimum hop routing. Note that the 
links in Fig. [T]are not only the Markov dependence 
edges but also the routing links. The minimum hop 
routing requires a hop count of |i — [n/2j | + |j — 
[n/2j I to deliver Yij to the fusion center. 

(A.4) The communication energy per link £'c((i„) ~ Egd'^, 
where > 2 is the propagation loss factor in wireless 
channel. 

(A. 5) Sensing requires energy, and the sensing energy per 
node is denoted by Es. Moreover, we assume that the 
measurement SNR increases linearly w.rt. Es, i.e., 
SNR = PEs for some constant (3. 

Henceforth, we consider various asymptotic scenarios and 
investigate the fundamental behavior of the ad hoc sensor 
network deployed over a correlated random field for statis- 
tical inference under assumptions (A.1)-(A.5). (Proofs are 
omitted due to limited space.) 

2 

The sensor density /i„ onX„ is given by pn = )2 ■ 

Assuming that the network is sufficiently large, the total 



information about the underlying field obtainable from the 
network is given by 

KLIt = n^Xs and MIt = n'^L, (24) 

and the total consumed energy in the network is given by 

n— 1 n— 1 

E = n^Es+EMY,Y.^\i-[nm\ + \j-[n/2\\), 
i;=o j=o 

= n^E, + Q{Tr')E,{dn). (25) 

Note that the knowledge of per-node information %s and 
Is and their properties w.r.t. SNR and sensor spacing d„ in 
(l24l i is critical for further development, and it is provided in 
the previous sections. 

We begin with the increasing area case. 

Theorem 3 (Infinite area and fixed density) For an ad hoc 

sensor network with a fixed and finite node density, the to- 
tal amount of information increases linearly as the area in- 
creases, but under both information measures the amount of 
harvested information per unit energy decays to zero with 
rate 

r) = e (area-^/^^ , (26) 

for any non-trivial diffusion rate a, i.e., Q < a < oo as we 
increase the area. 

Next, we consider the case in which the node density di- 
minishes, i.e., dn ^ CO. This case is of particular interest at 
high SNR since at high SNR less correlated samples yield 
larger per-node information. However, the per-sensor infor- 
mation is upper bounded as dn oo, and the asymptotic 
behavior is given by the following theorem. 

Tfieorem 4 As dn — > oo, the per-node information %s cind 
Is converge to 3C^(0) = D{Af{0, 1)||7V(0, 1 + SNR)) and 
/s(0) = ^ logSNR, respectively, and the convergence rate 
is given by 

K4d„) = -C4V^e-"'*" +o(v^e-"''") , (27) 

Is{dn) = ls{0) - c'^Vd^e'""'" + o (v^e-°'*") , (28) 

for constants C4, C4 > depending on the SNR. 

Theorem m can be proved using ( fTSl \\% and ( 1211 122] ). 
and explains how much gain is obtained from less corre- 
lated observations by increasing the sensor spacing in 2-D. 
Fig. |2] shows %s and E^ as functions of dn for a ~ \, 
€4 = 1 and 10 dB SNR. The gain in information is given by 
,y^g-Qd„ fQj. j^gg (^^^^ whereas the required per-link com- 
munication energy increases without bound, i.e., Ec{dn) = 
Eod'^ {v > 2). Since the exponential term is dominant in 
the gain as dn increases, the information gain obtained by 
increasing dn decreases almost exponentially, and there is 



no significant gain by increasing the sensor spacing further 
after some value. Hence, it is not effective in terms of en- 
ergy efficiency to deploy a very sparse network aiming at 
less correlated samples at high SNR. 




Fig. 2. Per-node information and per-link communication 
energy w.r.t. sensor spacing rf„ (SNR = 10 dB, a = 1, 

C4 = 1) 

The per-link communication energy can be made arbi- 
trarily small by decreasing the sensor spacing. To investi- 
gate the effect of diminishing communication energy E^. as 
dn 0, we now consider the asymptotic case in which 
the node density goes to infinity for a fixed coverage area. 
In this case, the per-node information decays to zero as 
d„ — > since ( ^ 1/4 as dn — > 0, and and Is{C) 

converge to zero as ^ ^ 1/4, as shown in Section [32] The 
asymptotic behavior in this case is given by the following 
theorem. 

Tlieorem 5 (Infinite density model) For the infinite den- 
sity model with a fixed coverage area, the per-node infor- 
mation decays to zero with rate 

3Cs^C5fin' +o{ti-'), (29) 

for some constant C5 as the node density fin — > 00. Hence, 
the amount of total information per unit area ( nats/m?) con- 
verges to the constant C5 as ^„ 00. Furthermore, in the 
case of no sensing energy, a non-zero energy efficiency rj is 
achievable if the propagation loss factor v = Z, and even an 
infinite energy efficiency is achievable if u > 3 as fin ^ 00 
for fixed area^ 

The finite total information for the infinite density and 
fixed area model follows our intuition. The maximum in- 

-Of course, this depends on the assumption of Ec{d„) = E(,d'^ for 
any dn > 0. However, this assumption may not be valid for small dn- 



formation provided by the samples from the continuous- 
index random field does not exceed the information between 
X{x, y) and Y{x, y) except for the case of spatially white 
fields. It is common that the propagation loss factor j/ > 3 
for near field propagation (i.e., dn 0). Hence, infinite en- 
ergy efficiency is achievable as we increases the node den- 
sity for a fixed area considering only communication en- 
ergy. Note that the total amount of information converges 
to a constant as we increases the node density. So, the infi- 
nite energy efficiency is achieved by diminishing communi- 
cation energy as dn 0. Considering the sensing energy, 
however, infinite energy efficiency is not feasible since we 
have in this case 



C5 + 0(l) 



i^Es + ein^-")' 



v>2, 
(30) 

as n — > CX3 for fixed coverage area. In this case the sens- 
ing energy rP'Es is the dominant factor for low energy effi- 
ciency, and the energy efficiency decreases to zero with rate 
O Thus, it is critical for a densely deployed sen- 

sor network to minimize the sensing energy or processing 
energy for each sensor. 

In the infinite density model, we have observed that en- 
ergy is an important factor in efficiency. Now we investi- 
gate the change of total information w.r.t. energy. We fix 
the node density and consider two scenarios to increase the 
required energy: One is to fix the coverage area also and 
increase the sensing energy, and the other is to fix the sens- 
ing energy and increase the coverage area. We assume that 
the network size is sufficiently large so that our asymptotic 
analysis is valid. The energy asymptotic behavior for two 
scenarios is summarized in the following theorem. 

Theorem 6 As we increase the total energy E consumed by 
a sensor network with a fixed node density and fixed area, 
the total information increases with rate 



Total information = O (log E) 



(31) 



as E oo. When the node density and sensing energy are 
fixed and the increasing energy is used to enlarge the cover- 
age area, on the other hand, the total amount of information 
increases with rate of 

Total information = 9 {e'^/^^ , (32) 

for any v > Q, as E oo. 

Theorem |6] suggests a guideline for investing the ex- 
cess energy. It is not efficient to invest energy to improve 
the quality of sensed samples from a limited area. This 
only provides the increase in total information in logarith- 
mic scale. Rather the energy should be spent to increase the 
number of samples by enlarging the coverage area even if it 
yields less accurate samples. 



5. CONCLUSIONS 

We have analyzed the asymptotic behavior of ad hoc sensor 
networks deployed over correlated random field for statisti- 
cal inference. Using our large deviations results that char- 
acterize the asymptotic information rate in 2-D for GMRFs 
under the CAR model, we have obtained fundamental scal- 
ing laws for total information and energy efficiency as the 
node density, coverage area and consumed energy change. 
The results provide guidelines for sensor network design for 
statistical inference about 2-D correlated random fields such 
as temperature, humidity, density of a gas on a certain area. 
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